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Abstract 

We propose a nonperturbative resummation scheme for the four-point connected quark-antiquark 

o 

Green's function G 4 that shows how the Bethe-Salpeter equation may be 'unquenched' with re- 
£>Y spect to quark-antiquark loops. This mechanism allows to dynamically account for hadronic meson 

decays and multiquark structures whilst respecting the underlying symmetries. An initial approx- 
imation to the four-point Schwinger-Dyson equation - suitable for phenomenological application 
is examined numerically in a couple of aspects. It is demonstrated that this approximation ex- 
plicitly maintains the correct asymptotic limits and contains the physical resonance structures in 



the near timelike region in the quark-antiquark channel whereas no resonances are found in the 
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diquark channel, respectively. 
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I. INTRODUCTION 



For many years, quantum chromo-dynamics (QCD) has been widely accepted as the 
underlying theory of the strong interaction. As an SU(N C ) gauge-field theory, QCD has as 
its elementary particle degrees of freedom the quarks and gluons, whereas the observables of 
the strong interaction are the spectra of hadrons and their spectral properties. One of the 
contemporary goals of QCD is clearly a detailed understanding of how the quark and gluon 
fields give rise to the formations of hadrons within a single theoretical framework. Whereas 
lattice calculations of QCD have shown substantial progress in the last years of dealing with 
dynamical fermions and improved actions [lj , final and robust results for light quark masses 
in the continuum limit are still lacking nowadays. 

In the absence of a final and complete solution to QCD one is motivated to adopt more 
pragmatic approaches to gain further insight to the strong interaction and its bound states. 
The Schwinger-Dyson equations are the field theoretical equivalent of the equations of mo- 
tion for the theory (see P, El 0,0] for recent reviews); they are dynamical equations in the 
continuum and embody all those symmetries that define the initial theory such as symmetry 
under Poincare transforms, the gauge symmetry, the discrete parity and charge conjugation 
symmetries and the (broken) chiral symmetry. The inclusion of all of these vital physical 
features gives an impetus to their study. However, the Schwinger-Dyson equations form 
an infinite hierarchy of coupled equations that must be truncated to some degree if actual 
numerical solutions are addressed. Contemporary truncation schemes have two guiding prin- 
ciples: i) the explicit maintaining of as many of the underlying symmetries as feasible and 
ii) agreement of the final results with observation (and possibly other approaches). 

The Schwinger-Dyson equations are coupled non-linear integral equations relating the 
Green's functions (as 'building blocks') of the theory to one another. By themselves these 
Green's functions do not have a physical interpretation but must be combined in various 
ways to construct physically observable quantities. One particularly efficacious framework 
jilQ, ^ is the Schwinger-Dyson equation for the quark propagator (gap equation) within 
the rainbow truncation as input into the ladder Bethe-Salpeter equation. The latter is itself 
a special case of a Schwinger-Dyson equation which leads to a description of the pion as an 
(almost massless) Goldstone boson - associated with chiral symmetry breaking - as well as a 
bound state of two (massive) constituent quarks. The mass of the constituent quarks here is 
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dynamically generated from the spontaneously broken chiral symmetry. The success of this 
construction in simultaneously describing two of the most fundamental aspects of hadron 
phenomenology can be traced back to the fact that the truncations employed explicitly 
observe the relationship imposed by the (flavor non-singlet) axialvector Ward-Takahashi 
identity (AXWTI) P, Q] . The AXWTI is an expression for gauge invariance when applied 
to quark-axialvector vector Green's functions. This shows the supreme role of symmetries 
when applied to dynamical systems. The light pseudoscalar and vector meson masses and 
leptonic decay constants 11, [3| as well as electromagnetic form factors Q are well 



reproduced with few parameters. In principle such parameters may be fixed by comparison 
with related quantities from lattice-QCD 15]. 

Whilst the success of the coupled rainbow-ladder Schwinger-Dyson - Bethe-Salpeter 
framework is laudable, it has so far proved difficult to go beyond this initial (simplest) 
truncation scheme. There have been several attempts to address different aspects of possible 
extensions. For example by using effective interactions characterized by 5— functions to 
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_ classes of dia™ U W and employing hybrid combinations of finite width and 
5— function interactions [18| . Ironically, it is the AXWTI that - preserving the symmetry - 
leads to a rapid increase of effort as one increases the level of sophistication in the truncation 
scheme. 

There are several issues of hadronic physics that the coupled Schwinger-Dyson - Bethe- 
Salpeter framework has not addressed so far. The first of these is the hadronic decay of 
mesons, e.g. the decay p — > titt. This question has been investigated within the impulse ap- 
proximation and acceptable results have been obtained for the decays of vector mesons ; 
however, the desired genuinely dynamic description is lacking. A second issue is the light 
scalar spectrum (cf. ref. (20 1 for an introduction to this topic). It is not clear presently, 
whether the lightest scalar mesons are simple quark- ant iquark resonances or if they are 
dominated by meson-meson, perhaps even diquark-diquark correlations. These issues are 
connected to unquenching, i.e. the inclusion of internal quark-loops to Bethe-Salpeter am- 
plitudes. The inclusion of such loops allows for multiple quark-antiquark or quark-quark 
correlations within the overall Bethe-Salpeter amplitude; these correlations give rise to 
the decay mechanisms and internal multiquark structures. Clearly then, such internal cor- 
relations play a significant role in the dynamical description of hadrons. The four-point 
connected quark-antiquark Green's function (G 4 ) is the simplest (and key quantity) of such 
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correlations. 

The aim of this paper is: i) to motivate the study of the 4-pt quark- ant iquark Green's 
function with a practical example for future more extended studies, ii) to show how a 
phenomenologically useful, dynamically generated approximation can be constructed, and 
Hi) to show actual numerical results. The example chosen is the unquenching of the Bethe- 
Salpeter kernel in a manner that preserves the AXWTI, maintains the charge conjugation 
properties and allows for the dynamical description of meson decay widths and possible 
multiquark systems. 

Throughout this paper we work in the isospin (d = u) limit with the only distinction 
between quark flavors being their current mass. We take two light flavors of quark, up and 
strange, and consider only the flavor non-singlet eigenstates uu, su and ss. This means that 
we do not consider the effects of (isosinglet and singlet) flavor mixing, such as the U(1)a 
anomaly induced in the r\ — r{ system, but rather focus on pure flavor eigenstates. We work 
in Euclidean space throughout, metric 5^ and with Hermitian Dirac matrices that obey 
{l^^lv} = 25 M „. The integral measure is dk = <i 4 fc/(27r) 4 . 

II. UNQUENCHING THE BETHE SALPETER KERNEL 

The 'unquenching' of a system entails allowing insertion of arbitrary numbers of internal 
quark loops into the various amplitudes of the system. Fully unquenching the system is 
clearly tantamount to solving a major component of the theory, which of course is beyond 
present techniques. However, as a first step we consider a restricted class of unquenching 
terms that include only a single quark loop and are diagrammatically planar. Hereafter, we 
refer to unquenching as the insertion of this single quark loop. 

As a prelude to unquenching the system we consider the fully amputated, connected 
quark- ant iquark (4-pt) Green's function G 4 that is two-particle irreducible with respect 
to any quark-antiquark pair under the ladder truncation. It obeys the Schwinger-Dyson 
equation (shown graphically in Fig. ^) 

G 4 Q/3;5 > + ,p_; k+) = Df u { P+ - M[^ a WK% 

+ J ^gG^. rK (p +) p_ ;g _,g + )[5 1 (g + )V; ] K7 [K 6 ^(?-)]^t(? + - K)- (1) 

In the above equation 
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FIG. 1: Schwinger-Dyson equation for the fully amputated, connected quark-antiquark (4-pt) 
Green's function, two-particle irreducible with respect to quark-antiquark pairs under the ladder 
truncation. Internal propagators here are fully dressed. 

Si{p) = -%MP 2 ) + s*(P 2 ) (2) 
denotes the quark propagator of flavor % and 

V; = iT a ^ (3) 

is the tree-level quark-gluon vertex. D a ^ v {q) and D a ^ u (q) are two different effective interac- 
tion terms, one for the seed term and one for the kernel of the equation. The reader might 
be worried about two different interactions terms, however, the reason for this separation 
will become clear in later sections and allow for practical truncation schemes that obey the 
underlying symmetries. Note, furthermore, that the distinction of seed and kernel interac- 
tions leads to an asymmetry between left and right quark-antiquark pairs. The momenta 
are given by p± — p ± Q/2, similarly for k and q. 

Equation (JIJ) - under the assumption that the G A contains resonant components in the 
timelike s-channel - can be used as a starting point for the derivation of the ladder truncation 
of the Bethe-Saipeter equation (see for example |2l|). As mentioned in the introduction 
earlier, these resonant components well describe the light pseudoscalar and vector mesons 
(the latter albeit in the case where meson decay is neglected). We emphasize though that 
Eq. ((H), unlike the Bethe-Salpeter equation, contains both the resonant and non-resonant 
components in a single expression. 

With the above 4-pt function (G 4 ) in mind we propose the truncated quark Schwinger- 
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Dyson equation 

Sai(p) =tt + m-fik [V«S(k)vXpDf u ( P - k) 

+ Jdk dqG^ Sl (p,p + q;k + q, k)[S(p + q)V«S(k + q)} p5 [S(k)V u % P Df u (P ~ k). 

(4) 

This equation is shown diagrammatically in Fig. 121 The first term of the self-energy is 
the standard rainbow truncation. The next term gives rise to unquenching since the 4- 
pt function G 4 (top line of Fig. |2J), when expanded (lower lines), generates planar graphs 
involving a single internal quark loop connected to the original quark line by successive 
numbers of gluon exchanges. By identification of the internal ladder graphs we collect those 
terms whose sum has been shown to give rise to physical resonances (at timelike momenta). 
The exterior interaction factors are given by D, whereas the purely interior interactions are 
provided by D; this will preserve the charge conjugation properties of the quark propagator 
(the left-right asymmetry of G A has been chosen accordingly). In the unexpanded form of 
the unquenching term, the explicit occurrence of the interaction D and the tree-level quark- 
gluon vertex V£ is to maintain the proper counting of the various graphs in the expansion. 
The explicit factoring of one tree-level vertex in each dressing term is a general feature of 
Schwinger-Dyson equations and is maintained here. Equations (jlj and JI}, once D and D 
have been specified, form a closed system which can be solved numerically and naturally 
includes possible resonance structures in the integrand at timelike momenta. One can regard 
the multiple gluon exchange terms in the unquenched quark Schwinger-Dyson equation as 
either a correction to the gluon propagator and/or vertex function or as a nonperturbative 
dressing of the quark by higher order correlations. 

The homogeneous quark-antiquark Bethe-Salpeter equation is 

r aP (p;P) = J dkK a ^(p,k;P) [Si(k+)T(k;P)S 2 (k-)] 7S (5) 

and the pole condition P 2 = — M 2 can be found once the kernel K and the quark propagators 
Si are specified. As emphasized earlier, the pion can be interpreted as the Goldstone boson of 
chiral symmetry breaking [l| when the kernel and the quark self-energy are related such that 
the AXWTI is satisfied. The flavor non-singlet AXWTI can be expressed in the following 
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FIG. 2: Proposed truncation of the quark Schwinger-Dyson equation to include a single quark 
loop with the use of the 4-pt function G 4 . All internal propagators are dressed. Springs represent 
the effective interaction D whereas wavy lines represent the interaction D (see text). 



There are two facets to this relation: on one hand, as a Ward-Takahashi type identity it is 
the expression of the underlying gauge symmetry written in a particular channel and is gen- 
uinely nonperturbative in nature; on the other hand it can be expanded semi-perturbatively 
to reveal relationships between different types of graphs. The former interpretation is re- 
sponsible for the chiral symmetry considerations and for the hope that - whilst one may not 
be solving the entire theory - one might be able to realistically describe a variety of physical 
phenomena. The latter interpretation allows to develop kernels K even though the relation 
introduced above is an integral relationship involving nonperturbative objects. 

We first introduce the kernel and show how it and the quark propagator from Eq. (j^j) 
satisfy the AXWTI, Eq. (JUJ). A short discussion of the physical implications will be given 
later on. The kernel is 



way [181 

[S-\p + )j 5 + l 5 S-\p^)] a(3 = [-i l5 f + 2m 75 ] Q/3 + J dkK(p. 



k ; P) a pfo [75 S ( k_ ) + S ( k+ ) 75] jS ■ 

(6) 




(7) 
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FIG. 3: Proposed truncation of the Bethe-Salpeter kernel to include a single quark loop. All 
internal propagators are dressed. Springs represent the effective interaction D whereas wavy lines 
represent the interaction D (see text). 

shown diagrammatically in Fig. El Just as for the quark self-energy this kernel contains 
the tree-level (ladder) term and a contribution involving G A . The expansion of G 4 into 
its diagrammatic form shows how a single quark loop is incorporated and each graph is 
planar. The kernel clearly has the correct charge conjugation properties. Knowing that 
the tree-level, rainbow quark and ladder Bethe-Salpeter terms already satisfy the AXWTI, 
Eq. (JHJ), we eliminate these terms. Using the diagrammatic expansion of the kernel and 
following the quark line, one sees that starting from left (index a) to right (index j3) or vice- 
versa on either side of Eq. (JHJ) the quark first and also finally interacts with the effective 
interaction D; all interactions between (if any) are with D. The quark loop, i.e. the only 
object proportional to the number of quark flavors Nf, is left untouched on both sides of the 
equation. This ordering and the fact that all planar diagrams are included imply that the 
AXWTI is naturally satisfied. This can be straightforwardly verified explicitly to arbitrary 
order without explicit specification of D or D. 
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Having established that the charge conjugation properties of and the AXWTI connecting 
the quark Schwinger-Dyson equation and Bethe-Salpeter kernel are satisfied, we now turn 
to the interpretation of the unquenching terms within the Bethe-Salpeter equation. There 
is the dual description of G 4 as: i) a sum of ladder diagrams and ii) as a quark- ant iquark 
correlation function. Recalling that the ladder approximated G 4 contains good represen- 
tations of the physical pions and kaons, the integral term of the kernel will now involve 
imaginary components picked up from integration over these resonances. This in turn gives 
the possibility for an imaginary component to the meson mass solution of the Bethe-Salpeter 
equation reflecting its finite hadronic decay width. In this way the dynamical description 
of meson decays may be achieved as well as a possible dynamical description of multiquark 
states and their mixing with the quark-antiquark components. The real contributions from 
the non-resonant parts of G 4 will only modify the position of the overall resonance. 

As a further remark, the form of the kernel makes clear that only flavor non-singlet 
combinations are included within this unquenching scenario. In fact, one can add more terms 
to the kernel that account for both the flavor singlet mixing and the anomalous breaking of 
chiral symmetry. Allowing for flavor singlet contributions one may 'open' the unquenching 
quark loop to form a kernel with two or more gluon exchanges in the s-channel. A similar 



system has been considered previously 
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23{ | . though with only two-gluon exchange. The 



kernel is still consistent with the AXWTI, up to anomalous terms. The reason for this 
is that while the quark self-energy can be constructed from functional derivatives of an 
effective action, the kernel is an additional functional derivative of the quark self-energy Q] . 
The flavor singlet contributions are then the functional derivative of the unquenching quark 
loop in our constructed quark self-energy Eq. (fl|). In this work, however, the kernel Eq. (J7J) 
only involves the functional derivatives of the original quark line and hence only involves 
flavor non-singlet parts. 

From the separate unquenching diagrams the quark loop is expressed as a trace over 
Dirac matrices; separately, each diagram is independent of the direction of the quark line. 
If the sum of these diagrams is replaced by the corresponding full 4-pt functions then again 
the direction of the internal quark lines is irrelevant. However, for approximate functions 
this is no longer true. We choose to use quark-antiquark 4-pt functions since we expect 
their resonances to form the basis for the physical decay of the overall meson. This outlines 
an important aspect of unquenching phenomena. We mention that with quark-antiquark 
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internal 4-pt functions one naturally expects an internal structure composed of mesons (and 
of course the original quarks) but the full 4-pt function contains also (non- resonant) diquark 
correlations. Moreover, one could in principle use quark-quark internal 4-pt functions and 
a diquark basis instead. The goal, therefore, is to maintain a reliable approximation to 
G A - in both resonant and non-resonant channels - such that the direction of the quark 
line does not affect results for physical observables. This picture of unquenching, where 
components are dynamically coupled to the quark-antiquark parts, serves as an illustration 
of the phenomenological descriptions of scalar mesons as qq 24J, meson-meson 25|, diquark- 
antidiquark [3, I^J or combinatoric [3] in character. In practice, all models can yield good 
results; the different descriptions essentially reproduce different parts of the full dynamical 
system. 

III. A TRACTABLE MODEL 

The Eqs. (JIJ), © and ((Zj) - once D and D have been specified - can be used in 
principle to find the resonance mass of a given channel while including a single internal 
quark loop. The major hindrance to solving this system is the number of components that 
must be kept track of. We here propose a suitable approximation that is expected to keep 
track of the most important pieces, and omit the rest for future study. 

The first step - in solving the unquenched Bethe-Salpeter system - is the specification 
of the effective interactions D and D. With the rainbow-ladder coupled Schwinger-Dyson - 
Bethe-Salpeter framework (which only uses D) , the Landau gauge and ultraviolet suppressed 
form 

DfM = 5 a %M)D(Q 2 ) = 5 a %M)^ 2 D^ exp (-^j (8) 

las proven to provide for a good description of the light pseudoscalar and vector mesons 
1^ | (in the absence of hadronic decays). The form of this effective interaction is tailored to 
primarily give results for hadron observables such as masses and weak decay constants. The 
success of the interaction Eq. (jHJ) can be attributed to the concept of integrated infrared 
strength - the explicit shape of the parameterization is unimportant compared to the integral 
area in the infrared region. One may add ultraviolet logarithmic terms to provide the correct 
perturbative limit at high momenta or compare with lattice calculations but these details, 
though important, are aot the 18 sue here (see e.g. Q for a more comprehensive discussion of 
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such issues). From the construction of the unquenched system it has become clear that D is 
the leading interaction term while D is sub-leading. Since the AXWTI can be satisfied even 
for a separate specification of the interaction terms, we have some flexibility in constructing 
the complimentary interaction D. Our initial criteria for choosing the form of D are that the 
resonance structure of G A for the lightest pseudoscalar mesons should be well approximated 
and that Eq. (0) is tractable numerically. We choose 



This form is similar to Eq. (jSJ) but is now in Feynman-like gauge. As an effective interaction 
it reproduces the leading structure of the pseudoscalar rainbow Bethe-Salpeter equation 
within Eq. ((TJ and retains the desired approximation to the tt and K mesons. It will also 
be seen that it allows for a Fierz reordering of the Dirac matrices of Eq. ; in this way the 
4-pt quark problem reduces to a single spin line. 

The next approximation is to neglect the back-reaction of the 4-pt function in the quark 
self-energy, i.e. we employ only the rainbow quark propagator. We stress that this is an 
initial approximation performed for heuristic purposes; the self-consistent set of Eqs. 
and (@J) will have to be examined in future in order to have a control on the approximations 
made here. In mitigation though, we posit that the major contribution of this term is to 
add possible imaginary parts and affect the resonance structure of the quark propagator 
functions in the timelike region. We expect that in the immediate vicinity of the spacelike 
axis (the region of interest here) there should be only small effects. 

The quark propagator is not itself an observable quantity and the precise details of its 
dressing are not directly of concern. One consequence of quark propagator dressing is the 
vacuum chiral quark condensate < qq >o7^ 0, which is a quantitative measure of dynamical 
chiral symmetry breaking |2.9]- The condensate < qq >o is calculable as an integral of the 
scalar part of the quark propagator at spacelike momenta alone and is real- valued. This gives 
credence to the assertion that the inclusion of the back-reaction term will not qualitatively 
alter the picture for spacelike momenta. Since < qq > is closely related to the properties 
of the pion, any minor quantitative changes could be compensated for by re-fitting model 
parameters. We also note that in [3fj] . where quark propagators constructed from simple 
complex conjugate poles have been investigated, the imaginary parts in the Bethe-Salpeter 
equation - generated by integration over these poles - canceled exactly, i.e. the kernel forbids 
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meson decay into quark constituents. Thus, although the details of the quark propagator in 
the timelike region will not be exactly correct when omitting the unquenched term, we do 
not expect the overall results to be seriously deficient. 
G 4 is decomposed into color components, 

where (for a given flavor combination input via the Si) G m and G n carry only Dirac and 
Lorentz structure. G m corresponds to a mesonic correlation whereas G n can be related to a 
diquark correlation. In the following we refer to G n as the diquark channel. 



Inserting the tree-level vertices, using the color identity 
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[ Ta W T 1<5/3 = 7^/3<^ 7 - Q S ^ S Sp, (11) 

the spin Fierz identity [32] 

brUhnhp = [K a U[K a ] S7 , K a = {1, * 7s , i lx /2, 2 7 57a/2} (12) 

and inserting the effective interaction forms Eq. (jSJ) and Eq. ® , the expression for the 4-pt 
function, Eq. (jTJ), becomes 



+i/ ^^((g-A:) 2 )G^(p + ,p_;g_,g + ) [5 1 (g + )^5 2 (g„)] r?K 



Ga/j;5 7 (p+,P-;fc-,fc+) = -3G^.*y(p+,P-; k-,k+) 

-5 / rfg^Gz - A:) 2 )G™ w (p +) p_; g_, g+) [5 1 (g + )Jf»5 a ( ff _)] IJ(( [K a 



57 
(13) 



where t = (p — k) 2 . Rewriting the solutions as 

G^(p +I p_; *_,*+) = G^(p + ,p_;A;_,fc + )[^], 7 -^(^-^) Q7 (^-^) 5/3 
G^> + ,p_;fc_,fc + ) = G^(p + ,p_;A;_,A ;+ )[^] 57 + ^(^-^) a7 (^-^)^ 



reduces the 4-pt equations to a problem involving a single spin-line 



(14) 



c££(p + ,p_;fc_,fc + ) = ^[r]^ 
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36 
1 



a q D((<i-kr D{(p - q)2) 



(p - 



a/3 



+ -J dqD((q-k) 2 )G n J(p + ,p^q_,q + )Tr[K b S 1 (q + )K a S 2 (q- 
3 J *^«*-*) a )-5^)T 



D{{ P - q y 



lit 



1 (q + )K a S 2 (q^(p l 



ia/3 



-3 J dqD((q - k) 2 )G^\p + ,p„;q.,q + )Tr [K b S 1 (q + )K a S 2 (q^)\ .(15) 

We stress that it is the choice of Feynman-like gauge for D that enables us to exploit the 
spin Fierz identity such that the Dirac structure of G m,n can be reduced with Eq. ()14j) . This 
reduction makes the analysis of G 4 feasible for phenomenological purposes. 

To proceed we make two further approximations (again for heuristic convenience). The 
first is to drop any terms in the trace proportional to the vectors q±\, i.e. we restrict to 
the leading covariant terms. This has the effect of suppressing any mixing between different 
channels, thus leaving only diagonal elements. In the spin-1 case only the terms proportional 
to the metric then survive. The second approximation is to drop the second term on the 
right-hand side of each equation. These last approximations allow us to write the original 
4-pt function and its Schwinger-Dyson equation in a much simplified form, which can be 
numerically evaluated (s = (p + — P-) 2 = Q 2 ): 



Gif}- Sl (p+,P-;k- } k + ) 



D(t) 



{$ - f) ai {i> - fig? 



+ [K a U[K a } 5j [5^5l y G m > a (s,t;k 2 ,k-Q) + 6^G n ' a (s,t;k 2 ,k-Q) 
where the (Lorentz scalar) functions G m,a and G n,a are the solutions of 
G n ' a (s,t;k 2 ,k-Q) = M + dqD((q-k) 2 )G n > a (s } t;q 2 ,q-Q)JC a 
G m ' a (s,t-k 2 ,k-Q) = -3G n ' a (s,t;k 2 ,k-Q)-^ I dqD((q - kf)G m ' a (s,t-q 2 ,q-Q)lC c 



(16) 



(17) 



with the kernels 



IC 1 
IC 5 

JC 5X 



4s 1 (q 2 + )s 2 (q 2 _) - 4g+-g_wi(g+)w 2 (g 2 ) 
-As 1 {q\)s 2 {q 2 _ ) - 4g + -g_w 1 (g+)t> 2 (g 2 ) 
si(ql)s 2 (q 2 _) - q+-q-Vi{q\)v 2 (q 2 _) 
Si{ql)s 2 (q 2 _) - q+-q~vi(ql)v 2 (q 2 _). 



(18) 
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There are several remarks in order here: The seed term in Eqs. (|17|) contains the only 
occurrence of the variable t. Since the solution will be proportional to this term, we can 
numerically set it equal to 1/ y/2 (see the next paragraph for an explanation of the Chebyshev 
convention) and thereafter express the solution in units of \/2D(t)/6. The function D is the 
original effective interaction dressing function and so the t dependence of G A is explicitly 
given by the single gluon exchange of the original seed term of Eq. (0) . We also note that as 
k 2 -> oo, G n ' a (s, t- k 2 , k-Q) -> D(t)/6, G m > a (s, t; k 2 , k-Q) -> -D(t)/2. This is seen by noting 
that D((q — A;) 2 ) vanishes exponentially as k 2 — > oo unless q 2 ~ k 2 — > oo in which case 
JC a vanishes. These considerations show that due to the form of the 4-pt dressing function 
the asymptotic limits of G A are preserved. Finally, the Bethe-Salpeter equation, albeit in 
approximated form, is still contained within Eqs. (|17|) with the position of the resonances 
dependent solely on s. This is important since we need to identify explicitly the l/(s + m 2 ) 
factors corresponding to meson or diquark propagation in order to specify the propagators 
for future applications. 

In order to solve Eqs. (fT7|) for a particular value of s, we perform a Chebyshev expansion 
in the angular variables z = k-Q/y/k 2 Q' 2 (and correspondingly z' = q-Q / \/q 2 Q 2 ) such that 

N-l 

Q(n,m),a^ f . ^ jfe.Q) (jH-«( Si £; fc 2 , Z ) = ^ ,m) '°(s, tj fc 2 ) T^z). (19) 

i=0 

Now the equations are a set of integral equations in a single variable {k 2 ). Note that we use 
a convention where T = l/v2. Equations (fTTj) are Fredholm equations of the 2nd kind, 
which can be solved for arbitrary complex s firstly for spacelike k 2 and then for general 
complex k 2 . 



IV. RESULTS FOR THE 4-PT FUNCTION 



In this section we present numerical results for the components of the 4-pt function 
G A . In what follows we use the parameters u> = 0.5GeV, D = 16GeV~ 2 , m u = 5MeV, 
m s = 115MeV (taken from 



12j ) and we set N = A. The (rainbow truncated) quark 



propagator functions are the same as in ref. [l2j. Looking at the form of the kernel Eq. (JJJ) we 
need to evaluate the functions G\ n ' m ^ a (s, t; I 2 ) at Q = P/2=fq (s = Q 2 ) and I = k+q/2±P/4, 
where P 2 = —M 2 is the complex mass solution of the Bethe-Salpeter equation, and k, q are 
spacelike momentum integration variables. 
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Feynman-like 


[12] 


J p 


uu 


su 


ss 


uu 


su 


ss 


0+ 


656 


957 


1230 


645 


903 


1113 


0" 


130 


446 


610 


137 


492 




1" 


912 


1540 
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TABLE I: Masses of mesonic resonances (MeV) in the 4-pt function G 4 . Results are compared 
with those of [12| which were calculated with all covariant terms, stable with respect to the number 
of Chebyshev moments (and shown to be explicitly Poincare covariant). 

The positions of the timelike singularities for the mesonic components G m (with the 
above parameters) are given in Table HI The (lightest) scalar and pseudoscalar results are 
surprisingly good, given the level of approximations applied. The vector and axialvector 
resonance positions are clearly deficient and this is especially severe in the asymmetric su 
flavor case. However, in our present study of unquenching effects in the light meson spectrum 
only the lightest resonances will be of importance due simply to the kinematics. 

The approximate rainbow-ladder Bethe-Salpeter meson masses (discussed above) are not 
the only singularities in G 4 . Also singularities for general complex values of s show up. In 
order to see the latter it is useful to consider the numerical solution of Eqs. (JT7J). When 
discretizing the momenta, these equations become matrix equations of the form T = C + KT 
and the condition for singularities is that Det(t — K) = 0, where K is the discretized 
integral kernel. A plot of \Det(t — K)\ as a function of complex s then shows the positions 
of such zeroes in each channel. The near timelike region studied is —lGeV 2 < dt(s) < 0, 
— lGeV 2 < ^s(s) < lGeV 2 while the spacelike region has been studied for 3?(s) — ► oo and 
increasing with dt(s). There are no resonances found in the spacelike region for any 
of the channels and also no resonances for the diquark channels, where \Det(t — K)\ is 
largely uniform and of 0(1). Figure 0] shows the zeroes in the near timelike region for 
the pseudoscalar and scalar meson channels with a uu flavor combination. In the near 
timelike region, the only other uu flavor resonance is the vector located on the negative real 
axis (as reported in the previous paragraph). The su flavor combination shows a similar 
behavior. It is not immediately clear what the extra complex conjugate pseudoscalar and 
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0,0 



0,0 



FIG. 4: Plots of |Det(l — K)\ for pseudoscalar (left panel) and scalar (right panel) channels as a 
function of complex s showing the location of zeroes in the near timelike region (uu case). 

scalar resonances are since they appear to be artifacts of the truncations and approximations 
employed. Anyhow, these resonances lie outside the region of interest for our studies of the 
light mesons. Nevertheless, the nature of these 'resonances' will have to be explored in 
future. 

Turning now to the solutions of Eqs. ()17j) we start with spacelike I 2 and show typical re- 
sults for the uu flavor combination with s = 1 + 1. As previously discussed, the t-dependence 
of the solution is removed by setting the seed term to 1/ y/2 such that the zeroth Chebyshev 
as I 2 — > oo tends to unity (—3) for the diquark (mesonic) components. We note that for 
general complex s the Chebyshev moments of the amplitudes will be complex. We show 
both the real and imaginary parts of the Chebyshev moments of the diquark G n ' 1 com- 
ponent in Fig. 03 Clearly seen is that the real part of the zeroth Chebyshev dominates 
primarily because of the seed terms in Eqs. ()17j) . Essentially the 4-pt function away from 
resonance positions shows its tree-level form. Furthermore, the diquark components do not 
have a resonance in either the spacelike or near timelike regions and one might conclude 
that G n could be well approximated by its tree-level term entirely; this will turn out not to 
be the case. The mesonic components for this value of s display the same behavior as their 
diquark counterparts. However, the mesonic sector does have resonances where the sub- 
leading Chebyshev structure and the deviation from the tree-level form become important. 
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FIG. 5: Real and imaginary (left and right panels respectively) parts of the Chebyshev moments, 



G^is = l + i,t,l 2 ) for real I 2 . 



We display the real parts of the zeroth Chebyshev amplitudes of the diquark and meson 
amplitudes as functions of spacelike I 2 for s = l + i in Fig. 13 This plot specifies the behavior 
of these functions away from resonance. As I 2 increases the scalar and pseudoscalar, vector 
and axialvector functions converge to the same values. This is because the difference between 
the channels (at this level of approximation) is the sign of the factor S1S2 in the kernel, which 
as I 2 — > oo rapidly vanishes. A comparison of the diquark and mesonic channels reveals that 
the two have similar shaped curves, though with different amplitudes. This can be attributed 
to the similarity in the kernels of the two channels. 

It is, furthermore, instructive to plot Go n,m ^ a (s, t, I 2 = 0) as a function of real s to highlight 
the resonance structure (cf. Fig. Cj). The diquark correlations show no resonance structure, 
but crucially do show significant variation with s, implying that the tree-level approximation 
is not reliable for low s. The mesonic correlations show resonance structures as expected. It 
is tempting to approximate these curves by a simple l/(s + m 2 ) shape, but closer inspection 
reveals that there they are indeed more complicated. This is an important point to note 
since the deviations from the pole term are significant. 

Finally, we show the functions Gf l ' rr ^ a (s, t, I 2 ) for the case of complex I 2 . Recall that 
in the kernel we need the amplitude for I 2 = (k + q/2 ± P/4) 2 , where k, q are integration 
variables, P 2 = —M 2 , while M is the complex mass of the meson under consideration. 
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FIG. 7: Zeroth Chebyshev moments of the diquark and meson (left and right panels, respectively) 

i(rj 



amplitudes G^ L,m ^ a (s,t,l 2 = 0) plotted as functions of real s. 



The factor P/4 does, however, mean that the deviation from the spacelike axis will only 
be slight. A consideration of the exponential interaction shows that the functions will be 
varying most significantly in the timelike and imaginary directions. Thus it is important 
that these deviations are taken into account, if only to assess their relevance to the final 
results. We show G™' 5 (s = 1 + t, I 2 ) (the leading Chebyshev moment of the tt) for complex 
I 2 in the vicinity of zero in Fig. |H1 It is seen is that the functions vary smoothly and - as 
expected - the variation becomes larger when progressing into the complex region. 
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FIG. 8: Real and imaginary (left and right panels, respectively) parts of the zeroth Chebyshev 
moments of the tt amplitude G™' 5 (s = 1 + «, t, I 2 ) plotted as functions of complex I 2 in the vicinity 
of I 2 = 0. 

V. SUMMARY AND CONCLUSIONS 

In this work, a proposed approximate form for the connected four-point quark- ant iquark 
Green's function (G 4 ) has been studied. The motivation for the study is an extension of 
the conventional Schwinger-Dyson - Bethe-Salpeter scheme to explore unquenching effects 
in the light meson sector. The expression for the approximated four-point quark- ant iquark 
Green's function has been derived and shown to be consistent with asymptotic limits. Fur- 
thermore, it reproduces the leading structure of the ladder Bethe-Salpeter equation for the 
(lightest) pseudoscalar mesons. The resulting equations have been solved numerically and 
typical results presented for both real and complex momenta. 

Most of the approximations employed in this work have been made for heuristic purposes 
and should be investigated in more detail in subsequent work. Nevertheless, our calculations 
demonstrate that the scheme proposed already gives acceptable results for the light meson 
sector and excludes resonances in the diquark channel. Clearly these initial results provide 
a basis for the study of dynamical meson decay and possible multiquark structures. 
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